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ABSTRACT 


This  Appendix  contains  complete  results  of  numerical  calculations  by  FE  code 
ABAQUS/Explicit  for  four  geometries  of  a  shear  layer  with  different  stress  concentrators. 
Those  four  geometries  are  shown  in  Fig.  1  of  this  Appendix. 

All  calculations  have  been  performed  for  VAR  4340  steel  with  hardness  ~50  HRC.  Several 
shear  velocities  have  been  assumed,  jfrom  10  m/s  to  150  m/s. 

Because  of  thermal  coupling  the  local  heating  around  a  stress  concentrator  leads  to  failure 
due  to  carrying  capacity  of  the  material.  Two  local  failure  criteria  have  been  introduced  into 
the  numerical  scheme.  The  first  one  is  based  on  the  maximum  inelastic  deformation  in  the 
adiabatic  conditions,  and  the  second  one  on  the  concept  of  the  strain  energy  density,  also  in 
the  adiabatic  case.  Both  criteria  give  similar  results  concerning  the  failure  patterns.  A  general 
conclusion  can  be  reached  that  during  impact  shearing  the  failure  is  more  sensitive  to  the 
adiabatic  stress  concentrators  at  higher  impact  velocities.  The  final  results  are  given  in  the 
form  of  graphs  and  tables. 
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1.  Geometry 


Four  different  specimen  geometries  were  studied  and  they  are  presented  in  Fig.  1.  The  loading 
conditions  are  similar  as  in  the  experiments  and  they  are  described  in  the  Final  Technical 
Report.  It  should  be  only  noted  that  there  exists  the  symmetry  in  the  numerical  problem,  so 
only  half  specimen  is  considered.  The  proposed  FE  models  are  demonstrated  in  Fig.  2  a,  b,  c 
and  d.  The  general  purpose  finite  element  code  ABAQUS/Explicit  was  chosen  for  the 
simulations. 


2.  Boundary  conditions 

All  boundary  conditions  are  shown  in  Fig.  2.  Part  A  is  loaded  with  a  constant  velocity  and 
part  B  is  supported.  Attachment  of  a  specimen  to  the  Hopkinson  tube  is  provided  by  the 
screws  that  do  not  assure  a  complete  fixation,  but  rather  serve  to  place  the  specimen  properly. 
Such  support  was  modeled  by  applying  contact  surfaces  between  the  specimen  and  the  tube. 
This  type  of  boundary  conditions  permitted  to  reflect  real  phenomena  observed  at  the 
interface  tube-specimen  like  small  local  plastification  of  the  specimen,  slight  rotation  of  the 
external  parts  during  impact  or  separation  between  the  specimen  and  the  tube.  The  tube  was 
modeled  as  an  infinite  elastic  half-space. 

Impact  velocity  V(t)  was  applied  at  the  upper  nodes  in  the  specific  way,  from  zero  to  the 
maximum  value  Vm-  The  total  rise-time  tjj]^  was  chosen  according  to  the  Gaussian  cumulative 

distribution  function  depending  on  the  impact  velocity  Vi .  The  mathematical  form  of  the  rise¬ 
time  function  is  represented  by 


F(x) 


=  1 


^min 


•exri 


(X-m)^ 

20^ 


Schematic  representation  of  the  rise-time  function  is  shown  in  Fig. 3.  The  constant  rise-time 
was  assumed  as  tjn=  2-10"^  s.  Values  of  parameters  m  and  o  in  are  m=0.5,  (7=0.125. 
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3.  Constitutive  relation 


A  constitutive  relation  was  sought  to  reflect  plastic  behavior  of  VAR  4340  steel,  ~50  HRC,  at 
high  strain  rates  and  at  wide  range  of  temperature.  The  explicit  form  of  the  constitutive 
relation  is  assumed  as  follows 


MT) 


B\ 


kToj 


(r.+rj^+T, 


^  T  V 
1— ^log-^ 

^  D  ^  T 


where  B,  po.  t),  n,  m  are  respectively,  the  modulus  of  plasticity,  the  shear  modulus  at 
T=300  K,  the  temperature  index,  the  strain  hardening  exponent  and  the  logarithmic  rate 
sensitivity.  To,  To  and  Tq  are  normalisation  constants,  and  p(T)  is  given  by 


p(T)  =  Po[1-A(AT)-C(AT)2];  AT  =  T  -  Tr 

where  x  and  T  are  respectively  the  shear  stress  and  shear  strain.  A,  C  are  constants. 

It  should  be  noted  that  application  of  the  constitutive  equation  is  restricted  to  T>300  K. 

The  strain  hardening  exponent  n  is  a  linearly  decreasing  function  of  temperature 

n(T)  =  no(l-T/rj 

where  no  is  the  strain  hardening  exponent  at  T=0  K  and  Tm  is  the  melting  temperature. 
All  constants  for  VAR  4340  steel  are  presented  in  Table  1. 


Table  1 .  Material  constants  for  VAR4340  steel,  ~50  HRC 


Const. 

Value 

Unit 

Const. 

Value 

Unit 

B 

1493.0 

MPa 

A 

5.047- 10-4 

1/K 

Ito 

80769.0 

MPa 

C 

1.036-10-7 

1/K2 

13 

0.39 

- 

D 

2914.0 

K 

n 

n(T)  in  text 

- 

-to 

622.0 

MPa 

no 

0.135 

- 

P 

0.9 

- 

m 

2.776 

- 

X 

38.0 

W/mK 

To 

300.0 

K 

P 

7890.0 

kg/m3 

To 

1.610-3 

- 

Cv 

460.0 

J/kgK 

fo 

106 

1/s 

Trr, 

1812.0 

K 
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A  failure  criterion  based  on  the  critical  inelastic  deformation  Fpc  was  introduced  in 
the  constitutive  relation.  Specific  value  of  Fpc  was  obtained  from  the  instability  point  3x^F=0 
on  the  adiabatic  t-F  curves.  A  failure  criterion  based  on  the  critical  energy  density,  where 
the  critical  plastic  work  dissipated  during  deformation  was  found  by  integration  of 
the  adiabatic  t-F  curve  up  to  the  instability  point  9t^F=0,  was  also  introduced.  Both 
approaches  give  similar  results  and  they  are  compared  in  Fig.  4. 

The  FE  simulations  were  carried  out  with  the  thermal  coupling.  Since  a  large  part  of 
the  plastic  work  is  converted  into  heat,  the  deformed  material  heats  up.  The  balance  of  energy 
with  the  heat  conduction  leads  to  the  following  relation 


aF 

at 


aT 


px^  =  pCv^-X 


a^T 

ay" 


where  y  is  the  axis  of  the  heat  conduction  and  X  is  the  heat  conductivity  (Fourier  constant). 
Material  constants  associated  with  the  process  of  thermal  coupling  are  also  given  in  Table  1. 

The  approach  to  the  thermal  coupling  which  is  formulated  above  permits  for  numerical 
analyses  of  all  temperature-coupled  problems  as  purely  adiabatic  as  well  as  with  the  heat 
conduction.  When  the  process  is  entirely  adiabatic,  ^  =  0  (no  heat  conduction),  the  heating  is 
uniform  in  an  elementary  volume.  Because  of  a  short  time  scale,  this  assumption  was  applied 
to  all  FE  simulations. 


4.  Simulation  data 

Several  impact  velocities  were  considered,  from  10  m/s  to  150  m/s.  Standard  geometry  was  of 
the  main  interest,  additional  impact  velocities  were  added  for  this  type  of  geometry  in  order  to 
observe  the  Critical  Impact  Velocity  in  shear  which  was  estimated  as  ~130  m/s. 

All  simulation  impact  velocities  are  given  in  Table  2. 


Table  2.  Simulation  velocities 


Impact  velocity 

Geometry  of  the  specimen 

Vi  [m/s] 

Standard 

U-geometry 

V-geometry 

Sharp  notch 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

80 

X 

X 

X 

X 

100 

X 

X 

X 

115 

X 

130 

X 

X 

X 

150 

X 
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5.  Tables  and  graphs  with  obtained  results 


The  values  presented  in  Table  3  were  found  from  the  graphs  presented  hereafter  to  give  a 
quick  comparison  of  the  values  obtained.  The  detailed  analysis  and  discussion  of  the  results  is 
given  in  the  Final  Technical  Report. 


Table  3.  Results  from  the  numerical  simulations 


Geometry 

V  impact 

n 

62 

r2 

Energy 

[m/s] 

[mm] 

[J] 

Q 

10 

49,4 

0,138 

0,069 

0,312 

0,156 

7,11 

30 

61,2 

0,181 

0,091 

0,362 

0,181 

6,68 

60 

62,4 

0,053 

0,027 

0,465 

0,233 

7,51 

80 

54,6 

0,037 

0,019 

0,322 

0,161 

5,05 

100 

55 

0,033 

0,017 

0,239 

0,120 

4,25 

115 

53,1 

0,035 

0,018 

0,241 

0,121 

48,3 

0,039 

0,020 

0,228 

0,114 

150 

28,7 

0,025 

0,013 

0,195 

0,098 

u 

10 

0,231 

0,116 

0,345 

0,173 

9,22 

30 

0,235 

0,118 

0,372 

0,186 

9,18 

60 

66,4 

0,304 

0,152 

0,381 

0,191 

10,69 

80 

67,1 

0,184 

0,092 

0,417 

0,209 

11,03 

100 

63,5 

0,119 

0,060 

0,366 

0,183 

8,61 

130 

52,3 

0,064 

0,032 

0,331 

0,166 

6,11 

V 

10 

47,9 

0,087 

0,091 

3,52 

30 

59,3 

0,111 

0,099 

3,92 

60 

52,6 

0,061 

0,031 

0,229 

0,115 

5,04 

80 

46,4 

0,046 

0,023 

0,247 

0,124 

5,43 

100 

44,9 

0,043 

0,022 

0,228 

0,114 

5,4 

130 

55,7 

0,059 

0,030 

0,241 

0,121 

5,29 

N 

10 

49,4 

0,053 

0,027 

0,079 

0,040 

1,37 

30 

55,4 

0,061 

0,031 

0,095 

0,048 

1,58 

60 

38,4 

0,021 

0,098 

0,049 

1,49 

80 

36,5 

0,011 

mumm 

0,091 

0,046 

Rmax  maximum  reaction  force, 

81  displacement  of  the  central  part  of  the  DS  geometry  at  the  beginning  of  failure  (elimination  of  the  first 
element), 

n  corresponding  nominal  shear  strain  (=6l/hs,  hs=2.0  mm), 

52  displacement  of  the  central  part  of  the  DS  geometry  at  the  final  failure  (the  last  element  eliminated). 
r2  corresponding  nominal  shear  strain  (=8l/hs,  hs=2.0  mm), 

Energy  dissipated  plastic  energy. 
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nOURE  CAPTIONS 


Fig.  1  Specimen  geometries  with  different  stress  concentrators. 

Fig.  2  FE  mesh  for  four  models;  a/  standard  geometry,  b/  U-geometry,  c/  V-geometry, 
d/  geometry  with  a  sharp  notch. 

Fig.  3  Graphical  representation  of  the  Gaussian  Cumulative  Function. 

Fig.  4  Criteria  of  failure  applied  in  FE  calculations. 

Fig.  5  Crack  propagation  in  standard  specimen;  Vi  =  10  m/s. 

Fig.  6  Crack  propagation  in  U  specimen;  Vj  =  10  m/s. 

Fig.  7  Crack  propagation  in  V  specimen;  Vi  =  10  m/s. 

Fig.  8  Crack  propagation  in  specimen  with  sharp  notch;  Vi  =  10  m/s. 

Fig.  9  Temperature  distribution  at  shear  layer  at  failure,  Vi  =  10  m/s,  a/  standard  geometry, 
b/  U  geometry,  c/  V  geometry,  d/  sharp  notch  geometry. 

Fig.  10  Plastic  equivalent  strain  distribution  at  failure,  Vi  =  10  m/s;  a/  standard  geometry, 
b/  U  geometry,  c/  V  geometry,  d/  sharp  notch  geometry. 

Fig.  1 1  Displacement  of  the  central  part  of  the  DS  geometry  at  the  beginning  of  failure 
(elimination  of  the  first  element). 

Fig.  12  Displacement  of  the  central  part  of  the  DS  geometry  at  the  final  failure  (the  last 
element  eliminated). 

Fig.  13  Dissipated  plastic  energy  vs.  impact  velocity;  4  types  of  specimen  geometry. 

Fig.  14  Reaction  force  vs.  time  for  the  chosen  impact  velocities;  4  types  of  specimen 
geometry;  a/  Vj  =  10  m/s,  b/  Vi  =  80  m/s,  d  Vi  =  130  m/s. 

Fig.  15  Maximum  reaction  force  as  a  function  of  impact  velocity;  4  types  of  specimen 
geometry;  a/  linear  scale,  b/  logarithmic  scale. 

Fig.  16  Nominal  shear  stress  vs.  nominal  shear  strain  curves  for  the  standard  geometry. 

Fig.  17  Nominal  shear  stress  vs.  nominal  shear  strain  curves  for  the  U-geometry. 

Fig.  18  Nominal  shear  stress  vs.  nominal  shear  strain  curves  for  the  V-geometry. 

Fig.  19  Nominal  shear  stress  vs.  nominal  shear  strain  curves  for  the  sharp  notch  geometry. 
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